TTP98-39 
PM/98-35 



hep-ph/9811275 



Short-distance tachyonic gluon mass and 1/Q 2 corrections 



K.G. Chetyrkinj], 

Institut fur Theoretische Tcilchcnphysik, Universitat Karlsruhe, 
Kaiserstrasse 12, 76128 Karlsruhe, Germany. 
E-mail : chet ©particle . physik . uni-karlsruhe . de 

| Stephan Narison 

Laboratoire de Physique Mathematique et Theorique 
O ^ UM2, Place Eugene Bataillon 

^ . 34095 Montpellier Cedex 05, France 

E-mail: narison@lpm.univ-montp2.fr 

m : 

V.I. Zakharov 

(N 

' Max-Planck Institut fur Physik, 

I/"") , Fohringer Ring 6, 80805 Munchen, Germany. 

' E-mail: xxz@mppmu.mpg.de 

00 
ON 

' Abstract 

We consider the assumption that a tachyonic gluon mass imitates short-distance nonperturbative 
CIh, physics of QCD. The phenomenological implications include modifications of the QCD sum rules 

for correlators of currents with various quantum numbers. The new 1/Q 2 terms allow to resolve 
in a natural way old puzzles in the pion and scalar-gluonium channels. They lead to a slight 
reduction of the values of the running light quark masses from the (pseudo)scalar sum rules and 
of a a (M T ) from r decay data. Analogously such terms only affect slightly the determinations of 
the running strange quark mass from e + e~ and r decay data. Further tests can be provided by 
precision measurements of the correlators on the lattice and by the e + e~ — » hadrons data. 
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1 Introduction 



In this paper we will outline phenomenology of nonperturbative short-distance physics in QCD based 
on introduction of a tachyonic gluon mass. 

To set up a theoretical framework, we will consider physical observables characterized by a large mass 
scale Q, Q 2 » Aq CD where Aqcd is the position of the Landau pole in the running coupling: 

a s (Q 2 )^ 7 , (1) 



(Q 2 /A 



/mIii [0 2 /A 2 QCL) 



where: 

bo = -^ = — (ll--n f ) (2) 

for rif flavours, is the first coefficient of the (3 function. Then in the limit of Q 2 — > oo, QCD reduces 
to the parton model while at finite Q 2 there are corrections of order (1/ lnQ) fc which are nothing else 
but perturbative expansions and power corrections, (Aqcd/Q) 71 which are nonperturbative in nature. 
To be more specific, concentrate on the correlators of currents J with various quantum numbers: 

n,/(Q 2 ) = i J e* x (0\T{J(x), J(0)}\0)dx, (3) 

where Q 2 = —q 2 and we suppress for the moment the Lorentz indices. Furthermore, to stick to the 
standard notations [0, || we introduce the Borel/ Laplace transformed II(M 2 ) where: 

Uj(M 2 ) ( -±Y n j{ Q 2 ) (4) 



(n- 1)! \dQ 

in the limit where both Q 2 and n tend to infinity so that their ratio M 2 = Q 2 /n remains finite. Then 
a standard representation for correlators (||) at large Q 2 goes back to the QCD sum rules |, | and 
somewhat schematically reads as: 



n. 7 (M 2 ) » (parton model) • ( 1 + ° J ^ + + O (In" 2 (M 2 /A 2 QCD ) , M" 6 ) ) (5) 

ln(M 2 /Ag CD 



where the constants aj, cj depend on the channel, i.e. on the quantum numbers of the current J(x), 
and we assumed that the currents have zero anomalous dimensions. Moreover the terms of order 
1/lnM 2 and M~ 4 are associated with the first perturbative correction and the gluon condensate, 
respectively. 

One of the central points about Eq. (||) is the absence of M~ 2 corrections. The reason is that Eq. (|) 
utilises the standard operator product expansion (OPE) and there are no gauge invariant operators 
of dimension d = 2 which could have vacuum-to-vacuum matrix elements. 

We will consider a modified phenomenology accounting for a correction due to a (postulated) non- 
vanishing gluon mass squared, A 2 . In this approach Eqs ^) are replaced by 



n,(M 2 ) a (parton model) • [ 1 + — ^- ^ + ^ + ^ + O (in" 2 (M 2 /A 2 QCD ) , M~ e '' 

QCD , 



ln ( M 2 /A 2 \ .U-' W 



(6) 

where bj are calculable in terms of A and depend on the channel considered. As we shall explain later, 
the introduction of the gluon mass is a heuristic way to account for possible existence of small-size 
strings. Moreover, it is easy to realize that in this case we expect a tachyonic gluon mass, A 2 < 0. Also, 
we have to confine ourselves to the cases when the A 2 correction is associated with short distances. 



2 Power corrections beyond the OPE 

The validity of the OPE beyond the perturbation theory cannot be proven without identifying more 
precisely nonperturbative mechanisms. One may discuss both infrared- and ultraviolet-sensitive non- 
perturbative contributions. As was emphasised rather recently |1 the general arguments based on 
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dispersion relations alone cannot rule out extra 1/Q 2 pieces. One can substantiate the point by the 
simple observation that the removal of the Landau pole from the running coupling: 



a ^q2^ 1 _^ 1 ^QCD ^ 

" U ' "6 ln(Q 2 /A^ CD ) ' fe ln(Q 2 /A 2 Qra ) ^ - A QCd) 



does introduce a Aq CD /Q 2 term at large Q 2 . Note that the ad hoc modification of the coupling in the 
IR region cannot be accommodated by the OPE. This simple example emphasises that the standard 
OPE is a dynamical framework based on a particular analysis of a generic Feynman graph, with large 
Q flowing in and out, and allowing some of the lines to become soft In this paper we accept 

the standard OPE to treat the infrared-sensitive contributions. 

As for the ultraviolet-sensitive contributions, an insight into power corrections is provided by renor- 
malons (for review and further references see 0). Renormalons are a particular set of perturbative 
graphs which result in a divergent perturbative series, J2 n ana s(Q 2 ) with a n ~ n\(— b^) n . Starting 
with n = N cr = l/boa s (Q 2 ) the product |a„a"| grows with n despite of the smallness of a s (Q 2 ). The 
most common assumption is that the rising branch of the perturbative expansion is summed up a la 
Borel: 

X, »!(-*>) a s (Q)^Jdt 1 + WQ2) . t qT- (8) 

Although some further arguments can be given to support this assumption ||, the Borel summation 
still looks an arbitrary solution and the theory is basically undefined on the level of Aq CD /Q 2 terms. 
Hence, one may argue that there should be nonperturbative terms of the same order, A.q CD /Q 2 which 
make the theory uniquely defined. Note that this kind of logic is very common in case of infrared 
renormalons (see, e.g., [S). There were a few attempts to develop a phenomenology of unconventional 
Aq CD /Q 2 corrections ||. Since there is no means to evaluate the nonperturbative terms in various 
channels directly the crucial ingredient here is the assumption on the nature of nonperturbative terms. 
In particular, a chiral-invariant quark mass and dominance of four-fermionic operators a la Nambu- 
Jona-Lasino have been tried for phenomenology Q . 

Here we will explore the phenomenology of the 1 /Q 2 corrections assuming that at short distances they 
can be effectively described by a tachyonic gluon mass. The assumption is of openly heuristic nature 
and can be motivated in the following way. Consider the static potential V(r) of heavy quarks and 
represent it as: 

V {r) « -^1 + kr. (9) 

Such a form of the potential is quite common at large distances r with k ~ 0.2 GeV 2 representing the 
string tension. At short distances one expects that the linear correction to the Coulomb-like potential 
is replaced by a r 2 term M. The vanishing of the linear correction at short distances is a reflection 
of the absence of the Q~ correction in the current correlators (see above). And vice versa, if the 
coupling a s {Q 2 ) has an unconventional Q~ 2 piece then the potential V(r) has a linear correction at 
short distances as well §|. 

The power-like behaviour of the gluon propagator was a subject of intense theoretical studies both by 
means of the lattice simulations (see, e.g., || and references therein) and by means of the Schwinger- 
Dyson equation (see, e.g., and references therein). However, these studies refer to large distances. 
As for the short distances, the smallest distances where the measurements of the potential on the 
lattice are available are of order O.lfm flTl|| . Remarkably, there is no sign of the vanishing of the 
kr term so far. Although it is not a proof of the presence of the linear term at r — ► since the 
measurements were not specifically dedicated to short distances, such a possibility cannot at all be 
ruled out and one is free to speculate that the linear term does not vanish at short distances. Thus, let 
us assume that the form ([)]) remains indeed true as r — > 0. There is still no obvious way to evaluate 
1/Q 2 corrections to other variables in terms of the same k . Note, however, that as far as the kr term 
is treated as a small correction at short distances, it can be reproduced by a tachyonic gluon mass A 2 
where @ 

- ^A 2 a k. (10) 
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In this framework A 2 is intended to represent the effect of small-size strings. From a calculational 
point of view, we simply replace the gluon propagator 

and check that the integral is saturated by large momenta, k 2 ~ Q 2 . One should notice that, to 
the approximation we are working, our analysis is gauge invariant. (An attempt to generalise the 
substitution in Eq. ( |IT|) to higher loops can be found in Ref. jyjj). 

There are obvious limitations to this approach. In particular, we do not include terms of order A 4 and 
do not go beyond one loop. Moreover, to use the Eq. (0) consistently at all distances one should 
have evaluated the anomalous dimension of the string tension k which goes beyond the scope of the 
present note (see, however, fl3|) . Also, the effect of strings at large distances cannot be reproduced 
by any modified propagator [|14| and this is one more reason to stick to short distances. Finally, the 
procedure described above assumes that at short distances the kr piece in the potential is associated 
with a vector exchange while at large distances it corresponds to a scalar exchange and there is no 
direct evidence to support such a change [|ll| . 

With these reservations in mind, it is clear that we would better concentrate first of all on the 
qualitative features of the approach with A 2 ^ 0. Let us emphasise therefore from the very beginning 
that there are qualitative effects around. Indeed, the estimate in Eq. ( |l0| ) implies a large numerical 
value for A 2 . Applying Eq. (|l0|), for example, at distances corresponding to a s = 0.3 we would find 
A 2 ps — 1 GeV 2 . Although this estimate cannot be trusted to the accuracy, say, better than factor of 
2 such a big gluon mass might seem to be easily ruled out since it would produce too much distortion 
on the standard phenomenology. We shall argue that it is actually not the case and A 2 » —0.5 GeV 2 
can well be admittable. 



3 Puzzles to be resolved 

Turning to the phenomenology, let us first note that not all the 1/Q 2 corrections in QCD are asso- 
ciated with short distances. Respectively, not all the terms proportional to A 2 can be used to mark 
short-distance contributions. For example, in case of DIS the 1/Q 2 corrections are coming from the 
IR region and perfectly consistent with the OPE. Thus, the class of theoretical objects for which an 
observation of the 1/Q 2 corrections would signify breaking of the OPE is limited. One example was 
the potential V(r) discussed above. Other examples are the correlator functions (|^) where the power 
corrections start with Q~ A terms. In terms of the expansion in A 2 , we are guaranteed then that terms 
proportional to A 2 are associated with short distances. If such terms arose from the infrared region, 
then the final result would contain nonanalytical in A 2 terms, A 2 In A 2 [|l5j and we shall see that it is 
not indeed the case. 

Even this limited set of variables exhibits actually a remarkable variety of scales which are not ex- 
plained by the standard OPE |Tq] . We will briefly review these puzzles since it is an important 
question, whether a novel phenomenology we are going to explore can resolve the puzzles of the ex- 
isting phenomenology. 

One of the basic quantities to be determined from the theory is at which scale the parton model pre- 
dictions for the correlators in Eq. (^) get violated considerably via the power corrections. To quantify 
the scales inherent to various correlators in Eq. (|j|) one may introduce the notion of M 2 rit which is 
defined as the value of M 2 at which the power corrections become 10 per cent from the unit term |16[ . 
While choosing just 10 per cent is a pure convention, the meaning of M 2 rit is that at lower M 2 the 
power corrections blow up. Moreover, the values of M 2 rit can be evaluated using either experimental 
or theoretical input. In the latter case the calculation is usually under control only as far as the power 
correction is relatively small and that is why M 2 rit is chosen to refer to a 10 per cent, not 100 per- 
cent correction. 

Consider first the best studied p-channel. On the theoretical side, H p is represented as: 

H p (M 2 ) « (parton model) ■ (l + + -^(^(G^) 2 ) + (12) 
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The value of M 2 rit is related then to the value of the gluon condensate, (^(G^) 2 ): 



M 2 rU {p - channel) r> y'lO^a.G 2 ) r> (0.6 - 0.8) GeV 2 (13) 

where the factor of 10 in the r.h.s. reflects the 10 per cent convention introduced above. An inde- 
pendent information on M 2 rit in this channel can be obtained by using the dispersion relations for 
n p (M 2 ) and integrating over the corresponding experimental cross section of e + e~ annihilation into 
hadrons, or by optimising the M 2 -dependence of the Borel/Laplace sum rules. The value of M 2 rit is 
then most sensitive to parameters of the p-meson and comes out about the same 0.6 GeV 2 as indicated 
by Eq. ([Hf). Since the gluon condensate parametrises contribution of the infrared region one may say 
that at least in the p-channel the breaking of the asymptotic freedom is due to infrared phenomena. 
Recent detailed experimental studies of the r-decays did not contradict this picture |l7]] . 
If one proceeds to other channels, in particular to the 7r-channel and to the ± -gluonium channels, 
nothing special happens to M 2 rit associated with the infrared-sensitive contribution parametrised by 
the gluon condensate. However, it was determined from independent arguments fl6|| that the actual 
values of M 2 rit do vary considerably in these channels: 

M c 2 rit (7r - channel) > 2 GeV 2 (14) 
M crit( 0± - gluonium channel) w 15 GeV 2 . (15) 

In the pion channel, the lower bound on M 2 rit can be determined by evaluating the pion contribution 
to the dispersive integral and equating it to the parton- model result. Since all other hadronic states 
give positive-definite contributions, at lower M 2 the parton model is certainly violated. In more detail, 
one gets: 



MU« - channel) > « 1.8 GeV 2 , (16) 



where the running light quark masses are evaluated at the scale M 2 , and f n = 93 MeV. A detailed 
M 2 optimisation analysis of the pion sum rule gives a value & : 

M 2 „ 4 (tt - channel) » (2 ~ 2.5) GeV 2 . (17) 

Although the difference of factor (3-4) in the values of M 2 rit in the p- and 7r-channels might seem not 
so big, it cannot be matched by the standard IR contributions. 

Note also that the difference in M 2 rit in the p- and 7r-channels was confirmed by the direct mea- 
surements of the corresponding correlators on the lattice | fL8| . These measurements also provide with 
an independent support to the idea that asymptotic freedom is violated at moderate M 2 ~ 1 GeV 2 
sharply, due to power corrections. It is most important in case of channels with large perturbative 
corrections, the fact which might give rise to doubts in relevance of the notion of M 2 rit 
In case of the scalar-gluonium channel, the value of M 2 rit follows |r| from a low-energy theorem for 
the correlator associated with the current a s (G 2 lv ) 2 . This low-energy theorem can be translated into 
correction to the parton model predictions at large M 2 : 

n G (M 2 ) « (parton model) ■ { 1 + ( f) ( -V + ... J , (18) 



boj \aj M 4 

where bo is defined in Eq. (^). Thus, the low-energy theorem brings in a large numerical factor 
12/(6o7r) (7r 2 /a 2 ) w 400 in front of the M~ 4 correction as compared with the p-channel. This changes 
dramatically the estimate of M 2 rit to: 

M 2 rit (0 + ~ gluonium) m 20M 2 rit (p - channel) m 15 GeV 2 , (19) 

and of the resonance properties, respectively (see Ref. [ [20, |l9|), in rough agreement with the in- 
stability analysis of the scalar |^0|, |l9) and pseudoscalar [2l|7l9) gluonia sum rules. Moreover, the 
huge power correction ( |l8| ) remains of a uniquely large scale and is not supported by any other large- 
scale correction which makes it very difficult to interpret the breaking of the asymptotic freedom in 
terms of resonances IfRl . 
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Also, the phenomenology build on the IR power corrections cannot resolve the problem of 77 -meson. 
As is known from general considerations [^2j a dynamical resolution of this problem is not possible 
without accounting for field configurations of nontrivial topology. On the other hand, the standard 
OPE is based on standard, i.e. perturbative Feynman graphs and cannot account for a nontrivial 
topology. Hence, one invokes direct instantons when the currents J interact at large Q 2 with small- 
size instantons |l6| . The instanton- induced contributions to the correlators in Eq. (||) were extensively 
studied in the literature (for review and references see J23|). We shall further comment on this point 
later. 

4 Some phenomenological implications 

Now, we add a new term proportional to A 2 to the theoretical side of IIj(M 2 ), see Eq.(^|). It is 
convenient to consider the effect of A 2 7^ channel by channel. 

4.1 Constraint on A 2 from the p-channel 

Phenomenologically, in the p-channel there are severe restrictions [p4[ on the new term b p : 



- (0.03 - 0.07) GeV 2 . 



(20) 



To find out the implications of this constraint for the gluon mass, we turn to the correlator of two 
vector currents Jy(x) = if>ij tl Tpj(x), viz. 



IMQ 2 ) = 1 J d 4 x e^^TJ^x) (J v ) f (0)|0>, 



(21) 



Here the indexes i, j correspond to quark flavours; m, is the mass of the quark i. To first order in a s 
and expanding in rriij we obtain: 



where 



and 



r(l) 



T (i) ,tt(1) 

l V,con > ll V,NO ' 



n 



V^con 



QQ 2 



(22) 
(23) 



(16^ 2 )H^ 



20 vn 
6- " 



3 ' " Q 2 Q 2 ^ Q 2 



m 

AI^q + 22—1 



9 9 

m m~L 



a s A 2 

Vq 2 



128 
T 



46 m 2 



rn 

Q 



2 

_ 2 

2 V<2 



-^+32C(3)- T - + 16 C (3)--18^ 



+ Q ^r l iQ - Q ^r l iQ + 6— l jQ - 



Q 2 



Q 2 



Q 



Q 2 



HQ 



(24) 



The above result is in MS scheme and the notations are as follows 

t 2 
12 



Li Tfl- 

m± =m,i rrij, I^q = log^), and l lQ = \og{-^). 



2 The results described in Eqs. [_' {[ ! l| |67| ) below have been obtained with the help of program packages MATAD 
Q and MINCER written in FORM 
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The normal-ordered quark condensates in Eq. ( p2| ) do not receive by definition any perturbative 
corrections and are displayed explictly only to make easier the discussion below. Note that the terms 
of order X 2 /Q 2 in Eq. (24) are \x independent and, thus, do not depend on the way how the overall 
UV subtraction of the correlator in Eq. (|2l|) has been fixed. Eq. (|24|) are written for the case 
of normal ordered quark condensate, that is the reason why quark mass logs appear there. As is 
well-known these are of long-distance nature and can (and even should) be absorbed into the quark 
condensate J2^, pOfl . To do this we should use normal non-ordered condensates, which means that 
perturbative theory contributions should not be automatically subtracted from the latter. Instead, 
these contributions are to be minimally renormalized. 

In order a s there are two diagrams (see Fig. 1) leading to nonzero vacuum expectation value of (ipitpi) 
in perturbation theory. A simple calculation gives 



■i r 



3m: 
to, A 2 



1 



4tt 2 



m(4: 

mi 



-5 + 61n- 



mi 



5 . (J, , 
3 to; 



(25) 



Now in order to find the polarization operator Ily ' in the case of normal non-ordered condensates one 



proceed by equating Hy\ to the combination Hy^oN 



Ily with the condensate values taken 



from Eq. 12a). The result of the treatment indeed does not contain any mass logs and reads 



(167r 2 )II 



V,NON 



20 m_ mi_ 



4Z 



7T 



f-16C(3) , " 7 '" 



4/ 



7T Q 2 
12 



o 

TO 



2 

2 Q 2 

2 



2 

ml 



24C(3)- 



16- 



m+ 



9 

m_ 2 

6 Q2-U 



128 , s 76 to , m 2 

-^ + 32 C (3)^ y - + 16 C (3)^ 



12— 



m~_ 

12 Q^ 



(26) 



(27) 



As it follows from Eq. (|25|) the quark condensate obeys the non-trivial RGE: 

2 d .—- . 3to„TOjA 2 a s 

where we have neglected all quartic quark mass terms. For further analysis, we shall choose the 
subtraction point fi equal to the sum rule scale M, as discussed in fl3l[ . 

In the light-quark case relevant to the p-channels we can neglect the m 2 terms and U p (M 2 ) 
simplifies greatly: 



n p (M 2 ) = (parton model) ■ (l + (^) 1 + (-^ + 8((3)) ) 



(28) 



or, in numerical form, 



iI p (Af 2 ) = (parton model) ■ I 1 



1 - 1.05 



M 2 



(29) 



This result was actually obtained earlier ]3^| . It is interesting that terms of order M~ 2 were looked 
for first in connection with proposal (see first paper in Ref. j6|) that UV renormalons can be 
imitated by a nonperturbative quark mass. In terms of the quark mass the correction is substantially 
larger and is —6m 2 /M 2 (see Eq. (p6|)). The conclusion of the analysis |24| was that the data allow 
only for an imaginary quark mass and quite a small one: 



-(71 - 114) 2 MeV 2 



(30) 
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(further analysis favoured rather lower end of the band ]33|). 

Now, that we are interpreting possible M~ 2 correction in terms of the gluon mass, the conclusion |2J] 
that a nonperturbative quark mass could only be imaginary fits very well our scheme with a tachyonic 
gluon mass. Moreover, the value of the quark mass, say, m 2 = —0.01 GeV 2 is translated into a much 
larger gluon mass, A 2 s» —0.5 GeV 2 . 

To summarize our discussion of the p-channel, wc find an important independent indication that the 
present data could be interpreted in terms of a tachyonic gluon mass: 



A 2 (l GeV) « -(0.2 - 0.5) GeV 2 , 



(31) 



where we have used (a a /7r)(l GeV) ~ 0.17 ± 0.02. This result indicates that, even relatively large 
masses , say, A 2 ps —0.5 GeV 2 cannot be ruled out. A further improvement of the (axial-) vector 
spectral function from e + e~ or/and r decay data would improve the present constraint on the eventual 
existence of such a term. 



4.2 Alternative estimate of A 2 from the 7r-channel 

We shall be concerned with the (pseudo)scalar two-point correlator: 

V> (5) (Q 2 ) = if d 4 x e^ x (0|TJ ( 6)(z)jf B) (0)|0), 

built from current of the bilinear light quark fields: 

J(5)(ac) = (m, + {-)m j )q i (i'y 5 )q i , 



(32) 



(33) 



having the quantum numbers of the pseudoscalar (O + ) it or scalar (0 ++ ) a°/<5 mesons. In the chiral 
limit (m u ~ m<j = 0), the QCD expression of the absorptive part of the correlator reads: 



1 3 

-ImV>( 5 )0) ^ (m f + (-) m .) —s 
7r an 



1 



a s \ ( „, ( s \ 17 A 2 



where one should notice that the coefficient of the A 2 term: 



(34) 



(35) 



In order to see the effect of this term, we work with the Borel/Laplace sum rule: 

f°° 1 
C(t) = J dte~ tT - lm^ (5) (t), 



and the corresponding ratio of moments: 



d 



n{r)=-—\ogC{r), 

where t< is the hadronic threshold. The QCD expression of the sum rule can be written as: 

3 



£( T ) = 77 m u (T) + m d (r) 



~ dr 
TZ(r) 



= 2 



8tt 2 
3 



1 + S + (5 2 + 5 x )t + 5 iT 2 + <5 6 t j 



8tt 2 

dC 

dr • 



m„(r) +Wd(r) Big] t 3 



'C, 



l + S , + l(5 2 +S x )r-^r 3 



where |34 



S = 4.821a s + 28.953a 2 
5' = 6.821a s + 57.026a 2 



(36) 



(37) 



(38) 



(39) 
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Table 1 : Ratio of moments TZ(t) and value of A 2 



tee l/M 2 [GeV- 2 ] C exp /{2flmi) R exp Rqcd -(a. A)A 2 [GeV] 2 



1.4 


1.20 


0.27 ±0.06 


0.66i° o j° 


09+ 08 

u - uy -0.06 


1.2 


1.29 


0.36 ±0.07 


0.79185 


o in +0 - 09 

u - lu -0.05 


1.0 


1.42 


0.46 ±0.10 




11+ - 07 

U.J.J-_o .06 


0.8 


1.61 


0.58 ±0.12 


1 10+ - 26 


0.12 ±0.06 


0.6 


1.89 


0.73 ±0.15 


1 24+ - 20 


0.12 ±0.06 


0.4 


2.29 


0.89 ±0.16 


i oq+0.15 
.13 


0.10 ±0.06 



while: 

S x = -4a s A 2 . (40) 
By keeping the leading linear and quadratic mass terms, one has E|: 

5 2 = -2(m 2 + m 2 - m„TO d )r, 

(5 4 = ^- j^(a s G 2 ) - (m d - ^y-) (™) + (u <-> , 

<5 6 = fc^7r 3 a s (uu) 2 . (41) 

TOi(r) and a s = a s (r)/7T are respectively the running quark masses and QCD coupling. The expression 
of the running coupling to two-loop accuracy can be parametrized as (y 2 = t" 1 ): 



a.(„) ee (^) = ai°) 1 - afgloglog^ + 0(a 2 ) , (42) 



with: 



a<°> ee 



1 



-Alog(^/A) 



(43) 



and (3i are the 0(a\) coefficients of the (3 function in the MS-scheme for rif flavours, which, for three 
flavours, read: 

01 = -9/2 , /3 2 = -8 . (44) 

A is a renormalization group invariant scale but is renormalization scheme dependent. The expression 
of the running quark mass in terms of the invariant mass rhi is || : 

™>i{v) = m l (-/3ia s (^))- 7l//3l |l±|(^-|)a s M+0(a 2 )| (45) 

where 7$ are the 0(a l s ) coefficients of the quark-mass anomalous dimension. For three flavours, we 
have: 

7i = 2 , 72 = 91/12 . (46) 
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We shall use in our numerical analysis, the QCD parameters compiled in 

A = (375 ± 50) MeV 
(a s G 2 ) = (0.07 ± 0.01) GeV 4 
(im) 1/3 (l GeV) = -(229 ±9) MeV 
k ~ 2-3 

A/ 2 = (0.8 ±0.1) GeV 2 (47) 



(uG^u) = M$(uu) . (48) 



where k indicates the deviation from the vacuum saturation assumption of the four-quark condensates, 
while M 2 parametrize the mixed quark-gluon condensate: 

K 

2 

One can notice that the ratio of moments 1Z(t) is not affected by the quark masses to leading order 
of chiral symmetry breaking terms. Therefore, we shall use it for extracting A 2 . We parametrize the 
pseudoscalar spectral function as usual by the pion pole plus the 37r contribution. The 3n threshold has 
been calculated using lowest order chiral perturbation theory, while two ir' resonances are introduced 
with the following widths and masses in units of MeV |37j : 

Mi = 1300 ±100 ri= 400± 200 and M 2 = 1770±30 r 2 = 310± 50 (49) 

Including finite width corrections, one obtains: 



ilnr0 5 (s) = 2/ 2 rr4 

7T 



1 S 

5(S- ml) + 0(s - 9?T 4) ( 167r 2J2)2 Y^Phad.(s) 



(50) 



where phad.{t) is given in j37j. We use, in our analysis, the largest range of predictions given by the 
two different parametrisations of the hadronic spectral function proposed in |37j, which correspond 
to the two values £ = -0.23 + i 0.65 (best duality with QCD) and £ = 0.23 + i 0.1 (best fit to 
the experimental curve for the observed 7r'(1770) in hadronic reactions) of the phenomenological 
interference complex parameter £ between the two tt' states. The result of this analysis corresponding 
to Tt exp (r) is summarized in Table 1, where we have also used t c ~ 2.5 GeV 2 as fixed from the duality 
analysis in |37). Therefore, we can deduce to a good approximation, at the r— stability region (r w 0.8 
GeV- 2 ): 

{a s /Tr)\ 2 ~ i [TZ exp - ftgjg] ~ -(0.12 ± 0.06) GeV 2 . (51) 

One can notice, by comparing this result with the one from the p-channel, that, to leading order in 
A 2 , the two estimates lead to the consistent value: 

(a s /7r)A 2 ~ -(0.06 - 0.07) GeV 2 A 2 (r w 0.8 GeV -2 ) ~ -(0.43 ± 0.09) GeV 2 . (52) 

One can also notice in Table 1 that the introduction of A has enlarged the region of QCD duality |37j] 
or t stability Q to a lower value of M 2 = 1/r of about M p . 

To summarize, the pseudoscalar channel can be considered as a success for the phenomenology with 
A 2 ? 0. 



4.3 Effects of A 2 on m u d and m s from the (pseudo) scalar channels 

Re-analysing the C sum rule by the inclusion of this new A term, one can also notice that the presence 
of such a term tends to decrease only slightly the prediction for the sum (rn u + rrid) of the running 
light quark masses by about: 

Sud PS -5.6% (53) 
Applying this effect to the recent estimates in |5?], [38|, on e obtains: 

(m u + m d )(l GeV) = (11.3 ±2.4) MeV => (m u +m d )(2 GeV) = (8.6 ± 2.1) MeV . (54) 
For the estimate of the strange quark mass from the kaon channel, one also obtains a decrease of: 

S us ~ -5% (55) 

which is slightly lower than the one of the pion channel due to the partial cancellation of the A effect 
by the m 2 one. 
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4.4 A 2 and the value of a s (M T ) from tau decays 

A natural extension of the analysis is to test the effect of A 2 on the value of a s extracted from tau- 
decays. For this purpose, we re-do the analysis from R T given for the first in [ ^9[ , by adding the new 
A 2 term. The modified expression of the tau hadronic width is (j^J : 

Rr = 3 ((\V ud \ 2 + \V US \ 2 ) S T EW [l + 5 T EW + S T PT + 5 T 2 +51 + ^ P }. (56) 

where: \V ud \ = 0.9751 ± 0.0006, \V US \ = 0.221 ± 0.003 are the weak mixing angles; S T EW = 1.0194 and 
S EW = 0.0010 are the electroweak corrections; 5% is a small correction from the running light quark 
masses; 8 T NP are non-perturbative effects due to operators of dimension > 4. The correction due to 

Table 2: Reduction of a s (M T ) for different values of A 2 
-A 2 [GeV 2 ] S^[%] 



0.16 


4± 1 


0.25 


7± 1 


0.34 


9± 1 


0.43 


11 ± 2 


0.52 


14 ±2 


0.61 


16 ±3 



the new dimension-2 term is: 

S > = 2 W (57) 

where, from Eq. (p9|) , b p = — 1.05a s A 2 and A 2 is given in Eq. (p2). At fixed order of perturbation 



theory, one has [|40|, |39|: 

S T PT = a s + 5.2023a 2 + 26.366^ + 0{a 4 s ) , (58) 

where we have truncated the series at the complete computed coefficients. We use the values of the 
QCD non-perturbative parameters given in Q and the experimental value |l7j R T = 3.484 ± 0.024. 
By comparing the extracted value of a s (M T ) obtained using A 2 = and A 2 ^ 0, we deduce in Table 
2, using the previous approximations, a reduction of the a s (M T )-value for different values of A 2 . The 
error given in Table 2 reflects the uncertainty on the eventual scale-dependence of the quantity a s A 2 or 
A 2 . The quoted central value comes from the average between the result corresponding to the freedom 
by taking one of these two quantities as the input number in the analysis. 



• If we use the value of A 2 as given in Eq. ([32|), then, we get to leading order in A 2 : 

<^~-(ll±3)%. (59) 

On the other hand, a comparison, at the r-mass, of the value of a s from the Z or/and the world 
averag e EJ with the one from r decays, where the latter includes corrections, gives a difference of 
about @7 

5%' ~-(5±10)%, (60) 

which is of the same sign and compatible with the correction in Eq. ([59|) . One should also have in 
mind that there are obvious limitations on the accuracy of the estimate in Eq. (|59|), due to the no 
control so far over terms proportional to A 2 and higher powers of a s (M 2 ), including the anomalous 
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dimension of the tachyon mass. These unknown effects might be taken into account by enlarging the 
error in Eq. ( |59| ) by about a factor 2. 



• Inversely, one can improve the previous determinations of A 2 from the p and ir channels, either by 
adding, in the different experimental fits jlTj] of the moments of the r-decays [[39| |4f|, the new A 2 term 
with the extra (compared with previous experimental fits) constraint that the sign of its contribution 
is unambiguously fixed, or by attributing the eventual small deviation of the values of a s from t and 
from e.g. Z— decays as being due to this new term. If one proceeds in the later way, and use the result 
in Eq. (|6^), one would obtain from Table 2: 

A 2 (Af T ) ~ -(0.2 ± 0.4) GeV 2 , (61) 

which is less conclusive than the result in Eq. (p2|). 



• Since the introduction of the gluon mass still presumably gives a better control over 1/M 2 terms 
which have been introduced so far ad hoc and are the major sources of theoretical uncertainties 
(4^]-|Q, one expects, within the present approach, an improvement of the theoretical errors in the 
determination of a s (M T ) by about a factor 2. 

• Furthermore, by comparing our previous results with some other approaches which resum the higher 
order terms of the perturbative series p3[-|i^|, the effect of the 1/M 2 term becomes manifest when 
one extends the analysis of the tau hadronic widths to a lower hypothetical r-lepton mass ||^, [l7| . 
We have therefore checked that until M T « 1.3 GeV, the agreement of the prediction with the recent 
data JlTj remains quite good. 

• Finally, one should notice that a naive redefinition of a s by a coupling which contains implicitly 
a A 2 /M 2 term, in the expression of the physical hadronic width R T , but not on the level of the 
propagator in Eq. (|Tl]), will lead to inconsistencies. 

4.5 Effects of A 2 on the value of m s from e + e~ and r decays 

In the determinations of m s based on the ratio |38| or on the difference [[l7| of the non-strange and 
strange quark channels, the leading terms which are flavour independent will not contribute. Taking 
the leading term in m 2 A 2 , which can be obtained from the previous expression of Ily , one can deduce 
that the presence of the A term implies a slight decrease of the value of m s of about (1-3)%, which 
is, however, smaller than the quoted errors (15%) in its determination. 

In the case of the determinations of m s from the alone AS — 1 part of the inclusive r decay rate 
one can check, after using the resummed series given there, and the expressions of the different 
corrections in the massless case given in the previous subsection, that the presence of A tends to 
increase by about (8-9)% the value of m s obtained from this method. Again, this effect is still smaller 
than the large uncertainties (about 30%) of this determination from this method. 
However, the effect of A tend to improve the agreement betwen the two different determinations. 

4.6 The quest of the scalar a°/S and iv^-channels 

Consideration of the correlator of the scalar Jg current (or of its strange analogue Jk*), where 

J s = —j= (uu - dd) (62) 
v2 

brings most interesting problems and challenges to the phenomenology with A 2 7^ 0. 

Indeed, in the limit of massless quarks, the (ordinary) QCD expressions of the pion and scalar sum 

rules only differ for the contributions of the dimension-six (or more) condensates, where here: 

^6 = -k^-^-n 3 a s (uu) 2 , (63) 

such that the strength of M 2 Ht or the optimization scale of about 2 GeV 2 [|| is obviously the same 
in both channels, i.e. much larger than the one of the p-channel. However, one should also notice 
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that, in this case of ordinary OPE (A 2 = and neglect of the instanton effects), the scalar sum rules 
can reproduce within a good accuracy the experimental mass of the a°/S meson, while the predicted 

In the same way, the estimate 



decay constant leads to a width consistent with the present data 

of the running light quark masses from the scalar sum rules [ [f9[ |2j, 35|, |36| within the ordinary OPE 
gives a prediction consistent with the analysis from other methods |4l| or from the e + e~ |58|, r decay 
data |Q or a global sum rule extraction of the light quark condensate . 
It is also obvious that M 2 rit associated with A 2 ^ is the same in the pion and <5-channcls: 

M 2 rit (-K — channel) ss M 2 rit (S — channel) . (64) 

What may be even a more drastic prediction is that the sign of the leading M~ 2 correction is the 
same in the two channels. Moreover, a detailed analysis of the ao/5 sum rule in the case A 2 ^ gives 
acceptable phenomenological results, and moves the scale of optimization to lower values of M 2 rit like 
the case of the pion. 

On the other hand, direct instantons predict the same M 2 rit in the both channels but the opposite 
signs for the leading correction (^3|. Thus, the scalar channel might be the right place to discuss in 
more details the choice between the direct instantons and the M~ 2 corrections. 

On pure theoretical grounds, both direct instantons and nonperturbative M~ 2 corrections seem to 
be indispensable parts of the QCD phenomenology. Indeed, instantons arc necessary to resolve the 
rj -problem. The nonperturbative M ~ 2 terms are necessary to match the M~ 2 uncertainty of the 
perturbative series due to the UV renormalons. 

All these qualitative arguments cannot fix, however, which correction (if any of the two) becomes 
important first at large M 2 . The signatures of the direct instantons and of the tachyonic mass are 
different in the 7T- and 5-channels. Namely, the tachyonic mass gives the same-sign correction while 
direct instantons result in opposite signs for the deviations from the asymptotic freedom. The existing 
lattice data Jl8[ indicate, to our mind, that it is rather a mixture of the two mechanisms which works. 
Indeed, the corrections in the two channels arc rather of opposite signs but the correction in the pion 
channel is substantially stronger. In the <5-channel the correction is much smaller than it should be if 
the instanton-model parameters are normalized to the 7r-channel data. 

Further data with better accuracy could be helpful. For example, it is not ruled out that in the 
(S-channel the correction is first positive because of the gluon mass and becomes negative at lower M 2 
because of the effect of the direct instantons. 

4.7 The gluonia channels 

We shall first be concerned with the two-point correlator: 

n G (Q 2 ) = ijd 4 x (0\TJ G (x) (J G (0)) f |0) (65) 
associated to the scalar gluonium current: 

j G = ^l(G^) 2 . (66) 



Its evaluation leads to: 



3A 2 

n G (M 2 ) = (parton model) ( 1 - + ... ) . (67) 



Thus, one can expect that the A 2 correction in this channel is relatively much larger since it is not 
proportional to an extra power of a s . With A 2 ~ —0.5 GeV 2 we obtain, by using the same 10 per 
cent convention as in section 3: 

M 2 rit (0 + - gluonium) w 15 GeV 2 (68) 
in amusing agreement with the independent estimate in Eq. (|l9|). 

Exactly the same phenomenon is observed for the case of the two point correlator of the pseudoscalar 
gluonium currents: the relative strength of the A 2 /A/ 2 term added to the parton result coincides with 
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that for the scalar gluonium displayed in Eq. (|67|). 



However, one should notice that a more quantitative analysis based on r-stability of the corresponding 
(pseudo)scalar sum rules |19]-|21 leads to a lower value of M 2 rit w (3 — 5) GeV 2 , but still much larger 



than the scale of the p meson. 

Let us consider now the case of the tensor gluonium with the correlator: 
€»pM ^*Jd 4 x e^(0|T^(:r)^(0)t|0) 

it ( q 2 



where 



+ ^2 ( -j9^9 P <y - qp.qu9 pa - q P q a 9pv + qp,q a g vp + q^q^g^p + 2^g P flw + q v q P g P c 

+ ~>P'o [g^gvp + g^gva-^g^gp^j , (69) 

Oftv = —G^Gya + —gp,vG a/3 G a P . (70) 

A direct calculation gives the following results for the structure functions ipf an d their respective 
Borel/Laplace transforms: 

^ " Tf^^^HIO' (71) 

V2 20 200 V 6 9 J 20 V 3 M 2 ) v ' 



^° 20 200 ^ V 4 12 J 20 V 2 M 2 ) y ' 

If, instead of considering 6^, we would introduce the total energy- momentum tensor of interacting 
quarks and gluons 6^ v , then various functions components of f/v^pc are related to each other because 
of the energy-momentum conservation. Indeed, requiring that 

tfwp<x% = 

we immediately obtain: 

^2=\Q 2 ^I and ^ = jQVl\ 
and, a consequence, the following representation of the function in Eq. (|69|): 

IplupM) = [VlipVua + VlurVvp ~ jfllifTlpaj ViQ 2 )-, (74) 



where: 

^ T (Q 2 ) = Q 4 -^I(Q 2 ), V^- 3 ^- < 7r »! 



To compare the sum rules based on A 2 ^ with the existing ones with A 2 = jljj |5lJ , one should have 
in mind that the results in Jl9| ] come from the sum rules corresponding to the correlator ip T defined 
in Eq. (f75j). Literally, these sum rules are not affected by the new A 2 terms because, in terms of the 



imaginary parts of the structure functions, the A 2 correction exhibited in Eq. (71) is proportional to 
\ 2 6(s) which vanishes once multiplied by an extra power of s. However, further analysis of once more 
subtracted sum rules as well as of the r— stability Jl9| would reveal now a large mass scale. 
To summarize, we conclude, from Eqs (|7l|-[75|) above, that the Borel/Laplace transforms of the func- 
tions ipf(i = 0, 2, 4) as well as tp T will have a scale analogue to the one of the (pseudo)scalar gluonium 
channels. This feature might reveal that all gluonia channels have an universal scale , in distinction 
from the instanton model |23) where the tensor-gluonium channel is expected to be insensitive to the 
largest scale exhibited in the spin-0 gluonium channel. Measurements of the corresponding correlators 
could provide therefore an interesting test of the theory 
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4.8 Heavy quarks: QQ and Qg-channels. 

Introduction of the gluon mass would result in a substantial change in the QCD phenomenology of 
heavy quark interactions. However, we expect that it would be rather a reshuffle of the parameters 
than a decisive test. Indeed, there arc successful description of the quarkonia states with the linear 
potential in Eq. (||) extrapolated down to r = 0, see, e.g., [j52|| . We plan to come back to these points 
in a future publication. 

5 Summary and conclusions 

To summarize our discussions of the various channels, the introduction of the tachyonic gluon mass 
allows to explain in a very simple and unified way the variety of scales of the violation of the asymp- 
totic freedom at moderate Q 2 in various channels. More specifically, the phenomenology based on the 
introduction of A 2 ~ —0.5 GeV 2 leads to the following successful predictions: 

• The correct sign and order of magnitude of the M~ 2 correction in the p-channel. 

• The new M~ 2 correction in the pion channel breaks the asymptotic freedom in this channel at the 
mass scale M 2 rit which is about factor of 4 higher than in the p-meson channel. This scale falls close 
to the value of M 2 rit determined independently from the values of / w and quark masses. The sign of 
the correction due to the tachyonic mass is also the one which is needed to bring QCD in agreement 
with phenomenology. 

• A natural explanation of M 2 rit in the scalar-gluonium channel (Eq. (pq)) which is much larger than 
M 2 rit in the p-channel and which was found independently from a low energy theorem (Eq. (|l9|)). 

• The account for the M~ 2 correction lowers by (11 ± 3)% the value of a s (M T ) determined from the 
r-decays, and may improve the accuracy of its determination by a factor of about 2. It also decreases 
by about 5% the value of the u, d and s running quark masses from the pseudoscalar channels. The 
presence of this term slightly decreases by about (1-3) % the value of m s obtained from e + e~ (3^] and 
t decay data p7| . It increases by about (8-9)% the results in the in (4f| from the individual AS = 1 
inclusive r decay channel. This effect, then, improves the agreement between the two different deter- 
minations. 

To overview the logic of our analysis, we were motivated to introduce a tachyonic gluon mass by the 
data on the QQ potential at short distances (see section 2). In this case A 2 ^ brings in a linear 
term kr at short distances which is, however, only a small correction to the Coulombic potential. Any 
precise measurement of this correction would require a subtraction of large perturbative contributions 
and as a result there are no error bars available on the value of k at short distances. Thus, the 
idea was to embed this (hypothetical) linear term in the potential into a relativistic framework and 
consider further consequences from this extension. Hence, the introduction of the tachyonic gluon 
mass. This extension allowed to analyse the current correlators. As a result, the bounds on the A 2 
narrow substantially. In particular, A 2 w — 1 GeV 2 which would naively correspond to the known 
value of k does not pass phenomenological hurdles, which is not so surprising. What is much more 
amusing, the value A 2 ~ —0.5 GeV 2 fits the data well (see above). 

Most remarkably, the introduction of A 2 brought large qualitative effects, namely a variety of the 
mass scales M 2 rit where the asymptotic freedom is broken at various channels. Various M 2 rit may 
differ by a factor of about 20. Further crucial checks of the phenomenology based on A 2 ^ could 
be provided by measurements of the correlators on the lattice. The et°/<5-scalar channel and gluonia 
channels appear most interesting from this point of view. 

However, the extension of the standard QCD sum rules by introducing A 2 ^ cannot resolve all 
the phenomenological problems. Namely, the rj -problem cannot be solved in this way and asks for 
introduction of instantons or instanton-like configurations. The existing data on lattice simulations 
might indicate that in the a°/<5-scalar channels, direct instantons become important at about the same 
M 2 as the effects of the gluon mass do. Further measurements of the current correlators on the lattice 
could provide checks of the phenomenology explored in the present paper. 
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